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We show that there is a way to unify distribution functions that describe simultaneously a signal in 
space and (spatial) frequency. Probably the most known of them is the Wigner distribution function. 
Here we show how to unify functions of the Cohen class, Rihacek's complex energy function, Husimi 
and Glauber-Sudarshan distribution functions. 

I. INTRODUCTION 

Distribution functions are widely used in optical physics (see [ll for a review) and in quantum mechanics where 
they are usually called quasiprobability distribution functions 0]). Probably the best known quasiprobability 

distribution is the Wigner function 0, H[ with applications in reconstruction of signals [f| and resolution @] in the 
classical world and reconstruction of quantum states of different systems such as ions [1] or quantized fields pi. M. Hoi. [Tll| 
in the quantum world. In this contribution we would like to re-introduce a lesser known quasiprobability distribution 
function, namely the Kirkwood-Rihaczek funtion [l], [l^, [H, Q EH > show how it can be related to the Wigner function, 
and express it as an expectation value in some eigenbasis, just as the other quasiprobability functions may be also 
expressed. 

II. BEST KNOWN QUASIPROBABILITY DISTRIBUTION FUNCTIONS 

A. Wigner function 

We start by introducing the Wig ner function, probably the best known. It may be written in two forms: series 
representation (see for instance [l6(), and integral representation 

W(q,p) = - J due""(q + -\p\q--), (1) 

for simplicity we use the Dirac notation here (see the appendix). In the above equation, p is the so-called density 
matrix. In 1932, Wigner introduced this function W(q,p), known now as his distribution function @, Q and contains 
complete information about the state of the system (p = \ip)(ip\). 

It may be written also as in terms of the (double) Fourier transform of the characteristic function 

w{a) = i / exp(a/? * _ ^W)^ ( 2 ) 

with a — (q + ip)/y/2 and where C((3) in terms of annihilation and creation operators is given by 

C(/3) = Tr{pexp(0 - (3*a)}, (3) 
also known as ambiguity function in classical optics (l7j . 

B. Q-function 

The Q or Husimi function [l8| , which is expressed as the coherent state expectation value of the density operator 

Q(a)= ^ / exp(a/T - a*/3)Tr{/5exp(-/3*a) exp(/3a f )}d 2 /3, (4) 



the alternative form is 



Q(a) = ~(a\p\a), 



(5) 
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C. Relating Q and Wigner functions in a differential form 

It is possible two group the Wigner and the Husimi functions 

F(a, s) = i J C(P, s) exp(a/3* - a*(3)d 2 p (6) 
where C(/3, s) is the characteristic function of order s 

C((3,s)=Tr{D({3)p}eMs\f3\ 2 /2) (7) 

with s a parameter that defines which is the function we are looking at. For s = 1 it is obtained the P-function, for 
s = the Wigner function, and for s = — 1 the Q- function. 
The Q-function is then 

Q(a) = J G(f3)exp(a(3* -a*(3)d 2 p (8) 

and for s = the Wigner function 

W(a) = J G(/3) exp(a/3* - a* f3)exp(\(3\ 2 /2)d 2 (3 (9) 

where 

G(J3) = ^Tr{P(/3)p} exp(- |/f /2), (10) 
The equation above may be written as an infinite (Taylor) series and inserted into ^ to obtain 

OO c\—Tl p 

W ^ = E — / G ^ ex p^* ~ a *® i^i 2n ( n ) 

n=0 ' ^ 

Considering the equality 



da da 

we can cast equation (jTTJ) into 



!> ° exp(a^* -a*/3) = - |/3| 2 exp(a/T -a*/?) (12) 



DO 



O — 77, / O Q \ W 

W = E^T 0(a), (13) 



or, finally [if 



n— v 



»'"')=«'-M»(-^^]^'«'- ll4 » 
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The analysis just done will help us to relate the Kirkwood-Rihaczck function with the Wigner function 
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D. Glauber-Sudarshan 

For the sake of completeness, we introduce another well-known distribution function: the Glauber-Sudarshan P- 
function (l9l [20|. First let us note that by using the coherent state eigenbasis, we can express the density matrix as 
the following double integral 

P=\ f I Hp\(3)\a)((3\d 2 ad 2 (3. (15) 



This representation involves off-diagonal elements (a\ p |/3) , and two inte gra tions in phase space. The next diagonal 
representation was introduced independently by Glauber and Sudarshan [l9L I20I ] 



p = j P{a) \a) (a\d 2 a (16) 

and involves only one integration. Using the equation above, we can write the Glauber-Sudarshan P-function in the 
form 

P(a)= J exp{a(3* - a*^)Tr{pexp(/3a t ) exp(-/3*a)}d 2 /3, (17) 



or 



P(a) = F(a, 1) = - / C(J3, 1) exp(a/T - a*(3)d 2 (3 (18) 
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E. Cohen-class distribution functions 

A function of the Cohen class is described by the general formula p], [22j 

W C = ^J J I <p(y + ^x')<p(y-^x')k(x,u,x\uy-^ vx '- u ' x+u ' v) (h:dx'd^ (19) 

and the choice of the kernel k(x, u, x' , v!) selects one particular function of the Cohen class. The Wigner function, for 
instance arises for k(x, u, x' , u') — 1, whereas the ambiguity function is obtained for k(x, u, x' , u')2n5{x — x')5(u — u'). 

III. LESSER KNOWN DISTRIBUTION FUNCTION: THE KIRKWOOD-RIHACEK 

QUASIDISTRIBUTION FUNCTION 

Now we turn our attention to a lesser known distribution, the Kirkwood-Rihaczek function, that may be written 
using the notation above as [24[ 

K((3) = J tPae? a *- fl * a e s ^r 1 C{a) . (20) 

This equation has been obtained from an equation similar to equations (0), Q and p7|) . i.e. taking the double Fourier 
transform 



K (q,p) = J dudve- mp e lvq Tr{pe iv ^e lu P} (21) 

and taking the trace as in equation (|56[) in the appendix. 

We will now do an analysis similar to the one done in subsection 2.D. We relate the Kirkwood-Rihaczek function 
to the Wigner function by using (|20p . via the following exponential of derivatives 

K(0) = e'^e^^W {13). (22) 
We now use the non-integral expression for the Wigner function (l6j 

W(0)=Tr\(-lfDHj3)pD(j3)], (23) 
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with D (/3) = e' 3 "* 13 a , the so-called Glauber displacement operator. We cast the above equation into the form 

W(P)=Tr[(-lfpD(20)], (24) 

where we have used the trace property Tr(AB) = Tr(BA) and the following identities(— 1)™ ((3) — D (/3) (—1)". 
Now we use the factorized form of the Glauber displacement operator [2l[ D (2/3) = e -2 '' 5 ' e 2/3at e~ 2/3 a to obtain 



2|/3|^ e 2/3aT e -2/Ta 



W(0) = Tr (-1) n pe 
Therefore, we have that the Kirkwood-Rihaczek function may be written as 



(25) 



Tr 



i o A i 



{-l) n pe 4 ^e 4 ^ T £>(2/3) 



The calculation of the exponential of derivatives of the Glauber operator will be tedious but straightforward, 
just write the main steps to obtain the final form, for instance, it is not difficult to show that 



By using that [21 



k=0 



we can express the above equation as 



2 00 ok 

e i^D {2(3) = H k (a f + a - (3*) ^-e s V 2/ra 



fc=0 



with Hk(x) the Hermite polynomials. From the above equation, is easy to obtain 

g2n 

r k=0 k=0 

and therefore 



- H k (x) tt = £ Hk + 2n 



— i 9 i 
e 4 g 



leWD (2/3) = £ £ ^f-^fc+a" ( St + S ~ ^) ^-e s V 2 ^ s 

n=0fc=0 



A"! 



Now we use the integral form of the Hermite polynomials *2 

oo 

2 P 



H p (x) 



{x + itf e~ f df 



(26) 
We will 

(27) 
(28) 

(29) 

(30) 
(31) 



(32) 



to obtain 



K(/3,(3*) 



e -/3* 2 e -2/3/3* 



OO oo 



dx / dte(- 2 ^ +2r+2/3 ) Jt e- 2a;2 e 2 ^( /3 * + «(x| e a2 e- 2/3 * a (-l) K p|x) (33) 



— oo — oo 



by using 



e~ lkt dt = 2ir6 (k) 



(34) 
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and taking k — 2y/2x — 2(3* — 2(3 we have 

oo 

K(/3,(3*) = 2^e-P' 2 e~ 2 W" J dxS (2y/2x - 2(3* - 2/?) e ~^ 2 ^{P'+p) (35) 

— OC 

x {x\e« 2 e~ 2fru {-lfp\x). (36) 
Making use of the identity 5 (ax) = we finally obtain 

KW,n=e-^\-^f^)'(^\e^(- i nf^) (37) 



or 



v2 



K (J3, (?) = X l\e^-^ (X\ Dt (-/?*) (-/?*) (-1)" p |X) (38) 



In this form, we have succeeded in obtaining the Kirkwood-Rihaczck function as an expectation value, in terms of 
position eigenstates, just as the Q-function, in terms of coherent states [see equation (5)], the Wigner and Glauber- 
Sudarshan functions in term of number states (l6j . 

IV. CONCLUSIONS 

We have shown that some distribution functions may be related through a method that allows the construction 
of some quasi-probability functions such as the Wigner, Glauber-Sudarshan and Husimi functions [23]. This method 
consists of obtaining the distribution functions from a double Fourier transform of an averaged exponential operator. 
If we use the exponential operator in terms of creation and annihilation operators we construct the already mentioned 
distribution functions. Via this method, but leaving the exponential operator in terms of position and momentum, 
and ordering (factorizing) the exponential in a convenient way, another function lesser used in classical optics, namely, 
Kirkwood-Rihacek's distribution function may be obtained. This function was recently introduced in quantum me- 
chanics by Praxmeyer and Wodkiewicz [3, [H| to have a phase representation of the Hydrogen atom. The connection 
between Glauber-Sudarshan and Husimi functions and functions of the Cohen class has been given, i.e. the adequate 
kernels. Finally, the Kirkwood-Rihacek has been given in term of an expectation value, in terms of position eigen- 
states, just as other distribution functions may also be given in terms of (sums of) expectation values. This may be 
of interest because it has been already exploited the fact that these forms allow reconstruction of quasiprobability 
distribution functions IIO, LLfll • 



V. APPENDIX 



In Dirac notation, we denote functions "/" by means of "kets" |/). For instance an eigenfunction of the harmonic 
oscillator [25| 

iMz) = ^L-^ffnCc). (39) 
V2 n n! 

is represented by the ket |n), with n = 0, 1, 2, ... In quantum mechanics, these states are called number or Fock states 
(see for instance Q). Any function can be expanded in terms of eigenfunctions of the harmonic oscillator 



f{x) = Y d C n xjj n (x) (40) 



n=0 

where 



f(x)Mx)dx (41) 
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and in the same way any ket may be expanded in terms of |n)'s 

oo 

l/) = E c » ( 42 ) 

n=0 

where the orthonormalization relation 

/oo 
lj>m( x )i>n(n)dx = $nm (43) 
-oo 

has been used. The quantity (m| is a so-called "bra". 

The basis set of kets \n) is a discrete one. However, there are also continuous basis. We can form one continuous 
basis for example with the function e lpq /V2n and the corresponding ket \p). First note that 

1 f°° 

W) = e- i{p - p )q dx = S(p - p% (44) 

so that 



or, in bra-ket notation we have 



rearranging terms we have 

\p) = (/ °°jp)mp)\p') = np'), (47) 

i.e. we have the completeness relation 

nOC 

\p)(p\dp=l. (48) 



e ipq /V2n= / 5(p-p'y p,l /V2ndp, (45) 



\p')= S(p-p')\p)dp= {p\p')\ P )dp, (46) 



Finally note that the function e lpq /V2tt is an eigenfunction of the operator — i-^ with eigenvalue p. For position, an 
"eigenket" of q is 

q\q) = q\q) (49) 

and an "eigenbra" 

ww = m- (50) 

We therefore find 

(gW-g)=o, (5i) 

that has as solution Q 

(q'\q) = S(q' - q). (52) 

We then can express the completeness relation 



1 = / \q)(q\dq, (53) 

J — OO 

such that 

\q)(q\\y)dq= / *{q)\q)dq (54) 

-oo J — oo 
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where ^(q) = (q\\^) = (q\^)- the density matrix p is defined simply as the ket-bra operator p = | ( M> | . 
The completeness relation serve us among other things to calculate averages, for instance 

/oo />oo 
\q)(q\dq\9) = / (Hf\A\q)(q\9)dq (55) 

or finally, 

/oo 
(q\*)(*\A\q)dq. (56) 
-oo 

Note that in the above equation we are simply adding "diagonal" elements, i.e. we have the trace of the the operator 
|^) (^| A. As the trace is independent of the basis, we can have it in terms of the discrete basis |n) 

oo 

= 5^<n|*)<*|i|n) = Tr{pA}. (57) 

n=0 
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